Low dimensional ordering and fluctuations in methanol-/3-hydroquinone-clathrate 

studied by X-ray and neutron diffraction 
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Methanol-/3-hydroquinone-clathrate has been established as a model system for dielectric order- 
ing and fluctuations and is conceptually close to magnetic spin systems. In X-ray and neutron 
diffraction experiments, we investigated the ordered structure, the one- dimensional (ID) and the 
three-dimensional (3D) critical scattering in the paraelectric phase, and the temperature depen- 
dence of the lattice constants. Our results can be explained by microscopic models of the methanol 
pseudospin in the hydroquinone cage network, in consistency with previous dielectric investigations. 



PACS numbers: 64.60, 61.10.Nz, 61.12.Ld, 77.22. Ch 

Reducing the dimensionality has a drastic influence 
on order-disorder phenomena, as it generally enhances 
fluctuations. Phase transitions are shifted to lower tem- 
peratures, or may be suppressed entirely. In the one- 
dimensional Ising model with nearest-neighbor interac- 
tion, thermal fluctuations suppress ordering at any fi- 
nite temperature, independent of the strength of the in- 
teraction. Theoretical models for low-dimensional or- 
dering and fluctuations are mostly tested on magneti- 
cally interacting materials, where the magnetic moment 
is well decoupled from the surrounding chemical lattice, 
and the interaction is short-ranged. Materials with low- 
dimensional, but long-ranged interactions like the dipole- 
dipole interaction are rare. Due to its small energy scale, 
the magnetic dipolc-dipole interaction is often negligible 
compared to exchange-type interactions, magneto-elastic 
interactions and anisotropies induced by spin-orbit cou- 
pling and crystal field. The electric dipole-dipole inter- 
action (EDD) is several orders of magnitude larger than 
its magnetic pendant. The electric dipole moment is con- 
nected to polar molecules or to displacements of groups of 
atoms against each other. Usually, the interaction with 
the surrounding lattice is far from being negligible. Mate- 
rials, which allow to test theoretical predictions for fluc- 
tuations and ordering of "pseudospins" interacting ex- 
clusively via the EDD interaction, are extremely rare. In 
this context, clathrates with polar guest molecules pro- 
vide a happy exception. They may undergo structural 
phase transitions driven by the interactions between the 
dipole moments of the guest molecules, with practically 
no coupling of the orientational degrees of freedom to the 
host lattice. 

The present article deals with methanol- f3- 
hydroquinone clathrate, where well defined methanol 
electric dipoles occupy the almost spherical hydro- 
quinone cavities. The electric dipoles arrange in chains 



'J ' *A N -C 




t hydrogen 
4 carbon 
• oxygen 



FIG. 1: The methanol molecule (CH3OH) in the quinone cage 
network. The atoms of the guest molecules are shown with 
their Van-der-Waals diameters whereas for the hudroquinone 
host lattice the stich-and-ball presentation is used. 



and are well decoupled from the surrounding host lattice. 
Thus the methanol clathrate can be understood in terms 
of interacting point-dipoles and has been established as 
a model system for dipolar ordering processesi*2*&^£. 

Figure ^ shows the methanol molecule, which carries 
a dipole moment of 1.69 DebyeS, in the hydroquinone 
cage. The shortest distance between the guest molecules 
is along the hexagonal c-axis. The geometry of the lat- 
tice combined with EDD interaction is responsible for 
the formation of ferroelectric chains along the c-axis, re- 
flected in the dielectric properties: Woll et al. have stud- 
ied the dielectric response as a function of temperature 
and frequencji^. The static dielectric constant along the 
hexagonal c-axis, e cs , shows a strong deviation from the 
Curie- Weiss behavior already at 250 K. The tempera- 
ture dependence of e cs at higher temperatures is well de- 
scribed by a quasi one-dimensional (ID) Ising model with 
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a ferroelectric coupling J c to nearest neighbors within 
the chain and a mean-field antifcrroclcctric inter-chain 
interaction with coupling constant Jji At lower T, 
e cs (T) deviates from this mean-field-Ising-behavior and 
finally decreases with decreasing T, indicating the on- 
set of three-dimensional (3D) antifcrroclcctric intcr-chain 
correlations. At T c =65 K the methanol-clathrate un- 
dergoes a transition into an antiferroelectrically ordered 
structure. 

While the dielectric experiments probe the macro- 
scopic response at Q — 0, at the center of the Brillouin 
zone, diffraction experiments give access to the corre- 
sponding underlying microscopic processes and allow to 
study the structural susceptibility at arbitrary points in 
reciprocal space. In this article, we first present X-ray 
diffraction results obtained in the ordered phase, which 
allow to develop the dipolar structure model. Then we 
turn to the paraelectric regime, where we studied the ID 
and 3D critical scattering by X-ray and neutron diffrac- 
tion experiments, as well as the temperature dependence 
of the lattice parameters. We compare our results to 
theoretical predictions and discuss them in the context 
of the preceding dielectric experiments. 



Experimental 



Single crystals of methanol-/3-hydroquinone clathrate 
have been grown from a saturated solution of quinol and 
methanol at 313 K. Completely deuterated samples were 
used for neutron experiments to minimize the incoher- 
ent background. Typical sample crystals were 2x2x2 
mm 3 . At room temperature, the methanol-clathrate 
has a rhombohedral R3 structure 8 . Throughout this 
article, we use the hexagonal notation, which is better 
adapted to the dielectric properties of the clathrate than 
the rhombohedral. Exceptions are marked explicitly by 
(. . ) r h- The hexagonal lattice parameters at 300 K are 
a=16.62 A, c=5.56 A. For the X-ray investigations, we 
used monochromatic CuK a radiation and a two-circle 
diffractometer, equipped with a closed cycle refrigera- 
tor. The single crystalline samples were oriented with c 
perpendicular to the scattering plane. This allows scans 
within the (hkO) plane. Scans within planes (hkl), I con- 
stant, I < 1 are possible by moving the detector out 
of the diffractometer plane. Rotating crystal exposures 
have also been taken. Here the crystal was rotated about 
c and the incoming beam was perpendicular to the c-axis. 
Diffuse ID critical scattering has been measured at the 
two-circle neutron diffractometer D23 (CEA-CRG) at the 
high flux reactor of the ILL in Grenoble, France. Here 
the ID axis, c, was mounted in the scattering plane. 




a)T=100 K 




b) T= 39 K 

FIG. 2: Rotating crystal photographs taken at T=100 K (a) 
and 39 K (b). The elliptical lines are due to scattering from 
the Al-windows of the cryo-refrigerator. 



Results 

Ordered phase 

Rotating-crystal X-ray photographs of methanol- 
clathrate were taken above and below the phase tran- 
sition (T c = 65 K). Figure |20 displays the photograph 
taken in the paraelectric phase at T=100 K, figure^) in 
the ordered phase, at 39 K. At high temperatures, hori- 
zontal layers of reflections (hkl) with / = 0, ±1,±2, . . . 
are visible. Below T c , additional reflections in layers 
with integer and half-integer values ofZ, Z = ±|,±|,... 
appear. The (hkl) values of the superlattice reflections 
can be determined from their position on the photo- 
graph. Careful analysis shows that modes with wave 
vectors k\ — (00^),^ and fc 2 = (\\\) r h condense at the 
phase transition on the boundary of the rhombohedral 
Brillouin zone at the Lifshitz points F (fci) and L (fe). 
To generate all superlattice reflections as linear combi- 
nations of k\ and ki, L and two arms of the F-star are 
required. This means, two of the rhombohedral lattice 
vectors double at T c . Additionally, the unit cell doubles 
in the hexagonal c direction (as indicated by the L-type 



3 



c 




FIG. 3: Projection of the 3D ordered antiferroelectric struc- 
ture onto the hexagonal basal plane: white circles: dipoles 
at z = 0, grey circles: z = |, black circles: z = — |. +/- 
indicates the orientation of the electric dipole. Below T c , two 
of the rhombohedral lattice vectors double (indicated by the 
dotted arrows). Along c, the dipoles order ferroelectrically. 
The ferroelectric chains arrange in ferro-ordered (100) layers 
with alternating sign from layer to layer. These layers are 
visualized by dotted lines. Solid lines indicate the hexagonal 
unit cells. 

vector k 2 = {\\\)rh)- 

The occurrence of k\ indicates the freezing of the 
methanol molecules into one of the six symmetry direc- 
tions of the R3 space group. The resulting structure of 
the dipole system in the dielectrically ordered phase is 
shown in Fig. [21 The dipoles are arranged in ferroelectri- 
cally ordered chains running along the hexagonal c-axis, 
which are arranged in sheets of alternating sign. As a 
property of the EDD interaction and the lattice symme- 
try the dipoles prefer an antiferroelectric arrangement 
perpendicular to the chains, which can not be satisfied 
for all chains, because of the triangular lattice symmetry. 
As a consequence, every dipole chain is surrounded by 4 
chains, which are oriented antiparallel, and by 2 parallel 
oriented dipole chains. The antiferroelectric interaction 
is thus satisfied by the majority. The 4:2 pattern fills the 
hexagonal basal plane in a regular fashion: Every dipole 
has the same environment of ferro- and antiferroelectri- 
cally oriented pseudospins. This 2-fci ground state is also 
obtained from Monte Carlo calculations^ and simple lat- 
tice sum calculations, considering the EDD interaction 
of the guest molecules, only (see Ref. 10] for the case of 
acetonitrile). 

The doubling of the unit cell along c, evidenced by 
the occurrence of L-type k 2 reflections, indicates an an- 
tiphase arrangement in neighboring cells along c. hi in- 
dicates a dimcrization of the dipoles in the ID chain, i.e. 
the appearance of pairs of strongly coupled pseudospins. 
Considering the EDD interaction in combination with a 
fiat potential in the cage center, the latter dimerization 
is energetically favorable, as we will discuss below. The 



ordering of the dipole system at T c is connected with a 
structural phase transition into the triclinic space group 
PI, connected with a small triclinic distortion of the host 
lattice below T^. 



Paraelectric regime 

The ID fluctuations in the dipole chains give rise to 
diffuse scattering which is distributed in planes perpen- 
dicular to the ID axis (and weighted by the methanol 
molecular form factor). The onset of 3D ordering close 
to the phase transition is indicated in diffuse scattering 
at the positions of the superlattice reflections (of the low- 
T phase). The correlation length along the dipole chains, 
£ c , and perpendicular, follows from the quasi ID Ising 
model, which has also been applied to fit the dielectric 
data, toi 

\c) \2J C ) 2|i|(V2) W 

and 

(ir) = (it) w™ (2) 

(t = (T — T c )/T c is the reduced critical temperature). 
Figure 0] shows diffuse ID and 3D scattering as measured 
with X-ray and neutron diffraction at different temper- 
atures in the paraelectric regime. We investigated the 
weak ID scattering on a neutron diffractometer because 
the contribution of the methanol molecules to the total 
scattering is higher in neutron scattering with deuterated 
samples than in the X-ray experiments. The correlation 
lengths £ c and £j_ follow from the width of Lorcntzian 
peak profiles and are shown in Figure |SJ The solid lines 
in the figure are fits after Eq. (JIJ and The fitting 
parameter J c is determined to J c — 202 K, in excellent 
agreement with the dielectric results (195 KA). The struc- 
tural data are well described by the quasi ID Ising model. 
The onset of the 3D correlations close to T c finally leads 
to the phase transition into a 3D antiferroelectrically or- 
dered structure. 

There is no indication of a static dimerization of the 
dipole chains in the paraelectric phase, i.e. no intensity at 
the superlattice positions of k 2 . But as a true ID feature, 
the dimerization could already appear in the dynamically 
ordered chain segments in the paraelectric regime, where 
the ID fluctuations already reach some hundreds of unit 
cells. We found no sign of diffuse scattering in planes 
with half- integer values of I, as caused by fluctuations 
of the dimerization above T c . From calculations of the 
methanol structure factoi* 1 ^ where we allowed a displace- 
ment of the dipoles along the ID axis of S = 0.25 A, the 
intensity of the diffuse scattering due to the dimerization 
can be estimated as being at least two orders of magni- 
tude smaller than that of the ID critical scattering. The 
corresponding correlation length £ sp should be smaller 
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FIG. 4: (a): Diffuse ID-scattering at Q=(3 4.5 1 + £) (neutron measurement) and (b) diffuse 3D-scattering at Q=(7 + £ 3 — 
0.43£ 0.5) (measured with X-rays). At 65 K the superlattice reflection of the ordered low-T phase is already visible as a sharp 
peak. (An Al-powder line has been cut around h = 6.55.) 
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T (K) T (K) 

FIG. 5: Correlation length £ c (a) und £j_ (b), parallel and 
perpendicular to the dipole chains as it was determined from 
Fig. 0] Solid lines are fits after Eq. Q and . Results from 
neutron diffraction are marked by (o), X-ray results by (□). 



than £ c , the corresponding critical scattering as a conse- 
quence broader. If the dimerization of the dipole chains 
already occurs in the paraelectric phase, this effect pre- 
sumably is too small to be seen in standard diffraction 
experiments. 

Figure shows the temperature dependence of the 
hexagonal lattice parameters a and c, as derived from the 
Bragg positions. Whilst the a-axis decreases linearly with 
temperature, c shows a distinct anomaly below about 
170 K. 




100 150 200 250 100 150 200 250 



T (K) T (K) 

FIG. 6: Temperature dependence of the hexagonal lattice pa- 
rameters a (a) und c (b) of the hexagonal cell. The dotted line 
in part (b) is a linear extrapolation of the high temperature 
linear behavior to lower T. 



Discussion 

The structure of the dipole system in the dielectrically 
ordered phase in methanol-/3-hydroquinone-clathrate, as 
determined from the analysis of the superlattice reflec- 
tions below T c , agrees well with the ground state calcu- 
lated by Monte Carlo simulations^ and lattice sum cal- 
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culations. The temperature behavior of the correlation 
length £ c and £j_ is well described by the correspond- 
ing quasi ID model by Scalapino et aii The struc- 
tural investigations therefore confirm the microscopic pic- 
ture for low dimensional ordering and fluctuations of the 
methanol pseudospins in the /3-hydroquinone cagework, 
as it was deduced from the dielectric experiments. 

An interesting point is the appearance of superlattice 
reflections of type "fc 2 ", because they indicate a doubling 
of the unit cell in the ID direction. The most proba- 
ble explanation of this doubling is a dimerization of the 
dipoles in the ID chains. Considering a ID chain of elec- 
tric dipoles on a deformable lattice (with distance d), in- 
teracting solely by the electric dipole-dipole interaction, 
a dimerization of the chain with respect to an alternat- 
ing distance d + S and d — 6 of the dipoles is energetically 
favorable. Because of the rather steep r -3 potential, a 
dipole can lower its energy much more by approaching its 
neighbor by a distance 8 than it at the same time looses 
potential energy by removing from the second neighbor in 
the ID chain. The small lattice distortion at T c (a first 
order structural transition) finally leads to a coupling 
of the two otherwise independent wave vectors k\ and 
&2. The freezing of the methanol- molecule and the static 
dimerization of the pseudospins along the chains there- 
fore fall together at T c . In a system with smaller coupling 
of the pseudospins to the surrounding host lattice, two 
successive transitions should occur. In ID S=| magnetic 
chain compounds, CuGe03 is a prominent exampleii, 
a structural phase transition into a statically dimerized 
structure is observed well above the magnetic ordering 
temperature. The nature of this spin-Peierls-transition 
is purely quantum mechanical. The pairing of neighbor- 
ing antiferromagnetic S = h spins to a S = singlet 
leads to an energy gap thereby suppressing thermal spin 
fluctuations and lowering the ground state energy. 

What is the reason for the unusual temperature be- 
havior of the ID lattice parameter c? Fig. [7| shows the 
anomalous part of the temperature dependence of the 
lattice parameter, i.e. the difference of the lattice param- 
eter c from Fig. b) to an extrapolated linear decrease 
with temperature and the measured and calculated cor- 
relation length £ c from Fig. [21 in the same plot. When 
choosing the right scaling factor, the anomalous part of 
c, Can, and the correlation length £ c fall on the same curve 
and both are well described by the quasi ID Ising model 
in Eq. (JTJ. That suggests a coupling of the ID fluctua- 
tions in the paraelectric regime to local strains. Such a 
coupling, however, is forbidden in the centrosymmetric 
space group R3; the coefficients of the (direct) piezoelec- 
tric effect all vanish^. The peculiarity of the present 
system is the unusual large amplitude of the paraelec- 
tric fluctuations. Because of the reduced dimensional- 
ity of the interactions, the ferroelectric fluctuations are 
strongly enhanced. The anomaly of c occurs below about 
170 K, where the correlation length already reaches about 
40 lattice units, £ c (170 K)/c ~ 40 (see Fig.EJ. Coupling 
via higher order terms in susceptibility and strain, which 
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FIG. 7: Temperature dependence of the anomalous part of 
the lattice constant c (•, see text for explanation) and the ID 
correlation length £ c (□). The solid line is the fit from Fig|S| 
after Eq. JTJ. 



become important due to the long ranged correlations, 
then presumably leads to the observed electrostriction. 
This effect is likely to occur in other low-dimensional 
systems, but has not been reported before. We do not 
know, why the electrostriction leads to a dilatation of the 
/3-hydroquinone host lattice (c an > 0) and a larger dis- 
tance of the pseudospins in the ID direction. Intuitively, 
a contraction along c should be energetically more favor- 
able. 

We investigated mcthanol-clathrates with methanol- 
occupancies less than 100 % to prove that the anomaly 
of c is a true property of the EDD interaction. By lower- 
ing the fractional occupancy of the cavities, the interac- 
tion between the pseudospins is diluted and the transition 
into a dielectrically ordered state shifts to lower tempera- 
tures. Samples with a larger percentage x of filled cavities 
(x > x c ~ 0.76) show conventional ordering via a first or- 
der phase transition into an antiferroelectrically ordered 
structure, whereas the others (x < x c ) freeze into dipolar 
glassesi*2i2i4. The temperature dependence of the ID lat- 
tice parameter c for several filling degrees down to 50 % 
is shown in Figure [S] Comparing the structural data 
with the static dielectric constant for different degrees 
of occupation, see e.g. Fig. 13 from Ref. 0, there is an 
analogous behavior: The anomaly is most pronounced for 
the completely filled methanol-clathrate. When reducing 
the filling, the anomaly shifts to lower temperatures and 
broadens. 



Conclusion 

The methanol clathrate can serve as a model system 
for dielectric ordering and fluctuations and is conceptu- 
ally close to magnetic spin systems. Because of the lattice 
symmetry and the peculiar property of the EDD interac- 
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FIG. 8: Hexagonal lattice parameter c for the concentrations 
x=0.97, x=0.79, x=0.73 and x=0.50. The curves are shifted 
vertically, for clarity. At 270 K, the lattice parameters for all 
concentrations coincide. 



tion, the latter clathrate has quasi ID dielectric proper- 
ties. From analysis of the superlattice reflections in the 
dielectrically ordered phase, we can confirm the theo- 
retically predicted 4:2 structure. The latter ground state 
can be calculated from Monte Carlo simulations and first 
principles, considering the EDD interactions between the 



guest molecules, only. We observe ID and 3D critical 
scattering in the paraelectric regime in neutron and X-ray 
diffraction experiments. The temperature dependence of 
the correlation lengths £ c and £ a is well described by the 
corresponding quasi ID Ising theory and agrees perfectly 
with results from dielectric investigations, thereby con- 
firming the microscopic picture of low dimensional or- 
dering and fluctuations of the methanol pseudospins in 
the hydroquinone cagework. The EDD interaction on a 
deformable lattice should be unstable against a dimcr- 
ization of the dipoles in the ID chains. At the dielectric 
ordering temperature T c we find superlattice reflections 
that indicate a doubling of the unit cell in the ID di- 
rection, i.e. a static dimerization of the dipole chains. 
We find no sign of diffuse scattering in half-integer lay- 
ers of I due to fluctuations of the dimerization. Below 
T=170 K an anomaly in the temperature dependence of 
the ID lattice parameter c occurs. The large amplitude 
of the ID fluctuations in the paraelectric regime leads 
to a coupling of the susceptibility to local strains and 
a piezoelectric effect in the otherwise centrosymmetric 
space group R3. This effect is likely to be observed in 
low dimensional interacting systems because fluctuations 
are strongly enhanced and long ranged. 
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